THE LAPLACE TRANSFORM OF THE CUT-AND-JOIN EQUATION AND 
THE BOUCHARD-MARINO CONJECTURE ON HURWITZ NUMBERS 

BERTRAND EYNARD, MOTOHICO MULASE, AND BRAD SAFNUK 

Abstract. We calculate the Laplace transform of the cut-and-join equation of Goulden, Jackson 
and Vakil. The result is a polynomial equation that has the topological structure identical to 
the Mirzakhani recursion formula for the Weil-Petersson volume of the moduli space of bordered 
hyperbolic surfaces. We find that the direct image of this Laplace transformed equation via 
the inverse of the Lambert W-function is the topological recursion formula for Hurwitz numbers 
conjectured by Bouchard and Marino using topological string theory. 
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1. Introduction 

The purpose of this paper is to give a proof of the Bouchard-Marino conjecture [3] on 
Hurwitz numbers using the Laplace transform of the celebrated cut-and-join equation of 
Goulden, Jackson, and Vakil [151 EZ|- The cut-and-join equation expresses the Hurwitz 
number of a given genus and profile (partition) in terms of those corresponding to profiles 
modified by either cutting a part into two pieces or joining two parts into one. This equation 
holds for an arbitrary partition ft. We calculate the Laplace transform of this equation with 
ft as the summation variable. The result is a polynomial equation |32j . 

A Hurwitz cover is a holomorphic mapping / : X — ► P 1 from a connected nonsingu- 
lar projective algebraic curve X of genus g to the projective line P 1 with only simple 
ramifications except for oo 6 P . Such a cover is further refined by specifying its pro- 
file, which is a partition ft = (fix > \xi > ■ ■ ■ > fie > 0) of the degree of the covering 
deg/ = \fj,\ = fix + • • • + fi£. The length i(fi) = i of this partition is the number of points 
in the inverse image / _1 (oo) = {px, ■ ■ ■ ,pe} of oo. Each part fi{ gives a local description 
of the map /, which is given by z i — > z^ in terms of a local coordinate z of X around 
Pi. The number h g ^ of topological types of Hurwitz covers of given genus g and profile fi, 
counted with the weight factor l/|Aut/|, is the Hurwitz number we shall deal with in this 
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paper. A remarkable formula due to Ekedahl, Lando, Shapiro and Vainshtein [61 [TBI EH] 
relates Hurwitz numbers and Gromov-Witten invariants. For genus g > and a partition 
\jl subject to the stability condition 2g — 2 + £(fi) > 0, the ELSV formula states that 

_ (2g-2 + %) + H)! ffijg /• AV(1) 

where is the Deligne-Mumford moduli stack of stable algebraic curves of genus g with 
t distinct marked points, (1) = 1 — ci(E) + • • • + (— l) ff c 5 (E) is the alternating sum of 
Chern classes of the Hodge bundle E on M g< e, ipi is the i-th tautological cotangent class, 
and Aut(/u) denotes the group of permutations of equal parts of the partition fi. The linear 
Hodge integrals are the rational numbers defined by 



<r m • • • r ntCj (E)) = _ Vi 1 • • • ^-(E), 

which is unless n\ + - ■ - + ne+j = 3g — 3+1. To present our main theorem, let us introduce 
a series of polynomials £ n (t) of degree 2n + 1 in t for n > by the recursion formula 

at 

with the initial condition £o(t) = t — 1. This differential operator appears in [17]. The 
Laplace transform of the cut-and-join equation gives the following formula. 

Theorem 1.1 ([32]). Linear Hodge integrals satisfy recursion relations given as a series of 
equations of symmetric polynomials in i variables t\, . . . ,ti: 

(1.2) Y, ( T ^giX))g^ ( ( 2 9 ~ 2 + i)L L (t L ) + J L l+ l{U)L\{i}{tL\{i})\ 
n L V i=l li J 

~ / , 2 / Vm T nL\{i,j} A g\ L )) ^-lW\{y}W\{y}J 



1 ' / 

i=l a,b \ 



€+1 



stable 



+ E {Ta-r r , J A^ i (l)) ai |/|+1 {T- f) T-„^A^ 2 (l)) a25|i7 | +1 j |a+l(ii)4+l(^)ln jLV{i} (*i\{i}), 
91+92=9 / 

where L = {1, 2 . . . ,£} is an index set, and for a subset I C L, we denote 

ti = (ti)iei, ni = {m\iel}, T n/ = ]^[r ni , = lj£ni(*0- 

T/ie last summation in the formula is taken over all partitions of g and decompositions of 
L into disjoint subsets I\J,J = L subject to the stability condition 2g\ — 1 + |/| > and 
2g 2 -l + \J\ >0. 

Remark 1.2. We note a similarity of the above formula and the Mirzakhani recursion 
formula for the Weil-Petersson volume of the moduli space of bordered hyperbolic surfaces 
of genus g with I closed geodesic boundaries [29, SOj • 

(1) There is no a priori reason for the Laplace transform to be a polynomial equation. 
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(2) The above formula is a topological recursion. For an algebraic curve of genus g > 
and I > 1 distinct marked points on it, the absolute value of the Euler characteristic 
of the ^-punctured Riemann surface, 2g — 2 + £, defines a complexity of the moduli 
space Mg,e- Eqn. (jl.2p gives an effective method of calculating the linear Hodge 
integrals of complexity n > from those with complexity n — 1. 

(3) When we restrict (|1.2|) to the homogeneous highest degree terms, the equation re- 
duces to the Witten-Kontsevich theorem of V'-class intersections [H [22j [55] . 

Let us explain the background of our work. Independent of the recent geometric and 
combinatorial works ESI EDI EZ] , a theory of topological recursions has been developed 
in the matrix model/random matrix theory community [TJ [Tl]. Its culmination is the 
topological recursion formula established in [llj . There are three ingredients in this theory: 
the Bergman kernel of an analytic curve C C C 2 in the xy-plane called a spectral curve, 
the standard holomorphic symplectic structure on C 2 , and the ramification behavior of 
the projection it : C — > C of the spectral curve to the x-axis. When C is hyperelliptic 
whose ramification points are all real, the topological recursion solves 1-Hermitian matrix 
models for the potential function that determines the spectral curve. It means that the 
formula recursively computes all n-point correlation functions of the resolvent of random 
Hermitian matrices of an arbitrary size. By choosing a particular spectral curve of genus 
0, the topological recursion [5J [Ti] [12] recovers the Virasoro constraint conditions for the 
V>-class intersection numbers (r ni • • ■ T ni ) due to Witten [38 J and Kontsevich [22j , and the 
mixed intersection numbers (k™ 1 /?™ 2 • • • r ni ■ ■ ■ T n< ) due to Mulase-Safnuk [31] and Liu-Xu 
[25] . Based on the work by Bouchard, Klemm, Marino and Pasquetti on remodeling the 
B- model topological string theory on the mirror curve of a toric Calabi-Yau 3- fold [2], 
Bouchard and Marino [3] conjecture that when one uses the Lambert curve 



as the spectral curve, the topological recursion formula of Eynard and Orantin should 
compute the generating functions 



of Hurwitz numbers for all g > and £ > 0. Here the sum in the first line is taken over all 
partitions n of length £, and Si is the symmetric group of £ letters. 



(1.3) 



C = {(x,y) | x = ye~y} C C* x C* 



(1.4) 






Figure 1.1. The Lambert curve C C C* x C* defined by x = ye 



-y 
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Our discovery of this paper is that the Laplace transform of the combinatorics, the cut- 
and-join equation in our case, explains the role of the Lambert curve, the ramification 
behavior of the projection tt : C —* C*, the Bergman kernel on C, and residue calculations 
that appear in the theory of topological recursion. As a consequence of this explanation, 
we obtain a proof of the Bouchard-Marino conjecture [3] . For this purpose, it is essential 
to use a different parametrization of the Lambert curve: 

x = e -(w+l) and y= 1 —±. 

y t 

The coordinate w is the parameter of the Laplace transformation, which changes a function 
in positive integers to a complex analytic function in w. Recall the Stirling expansion 

-kk k+n 1 ,„_! 
e —r, y= k n 2, 

which makes its Laplace transform a function of y/w. Note that the x-projection tt of the 
Lambert curve (|1.3p is locally a double-sheeted covering around its unique critical point 
(x = e _1 ,y = 1). Therefore, the Laplace transform of the ELSV formula (jl.lj) naturally 
lives on the Lambert curve C rather than on the u>-coordinate plane. Note that C is an 
analytic curve of genus and t is its global coordinate. The point at infinity t = oo is the 
ramification point of the projection tt. In terms of these coordinates, the Laplace transform 
of the ELSV formula becomes a polynomial in t-variables. 

The proof of the Bouchard-Marino conjecture is established as follows. A topological 
recursion of [11 j is always given as a residue formula of differential forms on the spectral 
curve. The Laplace-transformed cut-and-join equation (jl.2j) is an equation among primitives 
of differential forms. We first take the exterior differential of this equation. We then analyze 
the role of the residue calculation in the theory of topological recursion [3l [IT], and find 
that it is equivalent to evaluating the form at q £ C and its conjugate point q £ C with 
respect to the local Galois covering tt : C — > C near its critical point. This means all 
residue calculations are replaced by an algebraic operation of taking the direct image of the 
differential form via the projection tt. We find that the direct image of (jl.2p then becomes 
identical to the conjectured formula (jl.5p . 



Theorem 1.3 (The Bouchard-Marino Conjecture). The linear Hodge integrals satisfy ex- 
actly the same topological recursion formula discovered in [11] : 



(1.5) Yl ( T nTn L A V g (l)) g/+1 dL(t) ® di nL (t L ) 



n,n L 



A 9 ( l ))g,l P m{t, ti)dt ® dti ® d£ nLX{ . } (t L \ {i} ) 

j=l m,n L \ {l} 



+ ( E ( T « 7 * T »i A fl-l( 1 ))p-l^+2 



stable 



+ E E^ A 9i( 1 )) 91 ,i/i+i^ T ^ A ^( 1 )) 92 ,ui + i )p a , b (t)dt®di nL (t L ), 



I\JJ=L b ' n J 



where 



dL I {ti) = ®4rLMdU. 



iei 
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The functions P a ,b{t) and P n (t,ti) are defined by taking the polynomial part of the expres- 
sions 



P a , b (t)dt = - 



t 2ff_ i) (L+i(t)ib+i{s(t))+L+Mt))ib+i(t) 



P n (t,ti)dt ® dti = da 



ts{t) ( £ n+1 (t)ds(t) i n+1 (s(t))dt 



t - s{t) \ s(t) -U t-U 



where s(t) is the deck-tranformation of the projection n : C —* C* around its critical point 
oo. 

The relation between the cut-and-join formula, (jl,2p and (jl.5p is the following: 

Cut-and-Join Polynomial Equation Topological Recursion on 

Equation ~ , > on Primitives ( 11.2ft 7~i * Differential Forms HI. 51 1 
Laplace transform ' ' Direct image ' ' 



{Partitions} > Lambert Curve Galois Cover ^ 



Mathematics of the topological recursion and its geometric realization includes still many 
more mysteries [21 [Tl~l [27] . Among them is a relation to integrable systems such as the 
Kadomtsev-Petviashvili equations [11]. In recent years these equations have played an 
essential role in the study of Hurwitz numbers [20j [211 I2E1 EH ESI ES] . Since the aim of 
this paper is to give a proof of the Bouchard-Marino conjecture and to give a geometric 
interpretation of the topological recursion for the Hurwitz case, we do not address this 
relation here. Since we relate the nature of the topological recursion and combinatorics by 
a Laplace transform, it is reasonable to ask: what is the inverse Laplace transform of the 
topological recursion in general? This question relates Laplace transformation and mirror 
symmetry. These are interesting topics to be further explored. 

It is possible to prove the Bouchard-Marino formula without appealing to the cut-and-join 
equation. Indeed, a matrix integral expression of the generating function of Hurwitz num- 
bers has been recently discovered in [I], and its spectral curve is identified as the Lambert 
curve. As a consequence, the symplectic invariant theory of matrix models [TJ [9] is directly 
applicable to Hurwitz theory. The discovery of pQ is that the derivatives of the symplectic 
invariants of the Lambert curve give H g j(xi, . . . , X() of (11, 4h . The topological recursion 
formula of Bouchard and Marino then automatically follows. A deeper understanding of 
the interplay between these two totally different techniques is desirable. 

Although our statement is simple and clear, technical details are quite involved. We have 
decided to provide all key details in this paper, believing that some of the analysis may be 
useful for further study of more general topological recursions. This explains the length of 
the current paper in the sections dealing with complex analysis and Laplace transforms. 

The paper is organized as follows. We start with identifying the generating function 
(jl.4p as the Laplace transform of the ELSV formula (jl.ip in Section [21 We then calculate 
the Laplace transform of the cut-and-join equation in Section E] following [32], and prove 
Theorem 11.11 In Section [J] we give the statement of the Bouchard and Marino conjecture 
[3]. We calculate the residues appearing in the topological recursion formula in Section [5] for 
the case of Hurwitz generating functions. The topological recursion becomes the algebraic 
relation as presented in Theorem 1 1.31 In Section [6] we prove technical statements necessary 
for reducing (jl.2p to (jl.5p as a Galois average. The final Section [7J is devoted to proving 
the conjecture. 
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As an effective recursion, (jl.2p and (jl,5p caicuiate iinear Hodge integrals, and hence 
Hurwitz numbers through the ELSV formula. A computation is performed by Michael 
Reinhard, an undergraduate student of UC Berkeley. We reproduce some of his tables at 
the end of the paper. 
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In this section we calculate the Laplace transform of the ELSV formula as a function in 
partitions /j>. The result is a symmetric polynomial on the Lambert curve (11.31) . 

A Hurwitz cover is a smooth morphism / : X — > P 1 of a connected nonsingular projective 
algebraic curve X of genus g to P 1 that has only simple ramifications except for the point 
at infinity oo G P 1 . Let / _1 (oo) = {pi, . . . ,pf\. Then the induced morphism of the formal 
completion / : X Pi — > P^ is given by z i — > z^ with a positive integer //j in terms of a formal 
parameter z around pi 6 X. We rearrange integers /ij's so that \x = (fi\ > Hi > ■ ■ ■ > [i£ > 0) 
is a partition of deg / = = fix + ■ ■ ■ + m of length = t. We call / a Hurwitz cover of 
genus g and profile [i. A holomorphic automorphism of a Hurwitz cover is an automorphism 
<fi of X that preserves /: 



Two Hurwitz covers /i : X\ — > P and /2 : X2 — > P are topologically equivalent if there is 
a homeomorphism h : X\ — » X2 such that 



2. The Laplace transform of the ELSV formula 



X 



X 




It 



x 2 




pi 



The Hurwitz number of type (g, fj) is defined by 
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where the sum is taken over all topologically equivalent classes of Hurwitz covers of a given 
genus g and profile fi. Although h g « appears to be a rational number, it is indeed an 
integer for most of the cases because / has usually no non-trivial automorphisms. The 
celebrated ELSV formula [6j [TBI [35] relates Hurwitz numbers and linear Hodge integrals 
on the Deligne-Mumford moduli stack M. g ,t consisting of stable algebraic curves of genus 
g with £ distinct nonsingular marked points. Denote by ir g ^ : Mg^+i — * M g ,e the natural 
projection and by uj n e the relative dualizing sheaf of the universal curve ir gi £. The Hodge 
bundle E on M. g ,i is defined by E = {Fgj)*tjJn e , and the A-classes are the Chern classes of 
the Hodge bundle: 

A; = Cj(E) e H 2i (M g/ ,Q). 

Let Oi : Mg t e — ► M- g ,t+i be the i-th tautological section of n, and put £j = a* (u^ t ) . The 
V'-classes are defined by 

^ = C1 (A) g H 2 (M g>e ,Q). 

The ELSV formula then reads 

, _ H fiVf y A 3 V (1) 

- iAuto*)i fi^^ (M) nS(i-^)' 

where r = 2g — 2 + + \fi\ is the number of simple ramification points of /. 

The Deligne-Mumford stack Ai g £ is defined as the moduli space of stable curves satisfying 
the stability condition 2—2g—£ < 0. However, Hurwitz numbers are well defined for unstable 
geometries (g,£) = (0, 1) and (0, 2). It is an elementary exercise to show that 

U _ t fc-3 flTlH u _ (Ml +M2)! Ml* 

fto.fc - « and n>o,(ui,u 2 ) - ; — r P 

H1 + H2 Ml' M2! 

The ELSV formula remains true for unstable cases by defining 

A V (1) . . 1 



A V (1) 1 



(2.1) 
(2.2) 

/ATo,2 i 1 ~ Mi^iX 1 - ^2) Ml + M2 
Now fix an £ > 1 , and consider a partition /i of length £ as an ^-dimensional vector 

n = (fi u . . . , fi e ) eN £ 
consisting with positive integers. We define 

|AutQz)| 
(2g-2 + e + n 1 + -- + /i t )\ 

A^(l) 



(2-3) g g (M) = , n . n /r r: — — 



Mi' 



^ fH P AVI 

= n 1 n' ■ L = £ <T - • • • r "' A » v<1>> n - 

CIS £1 function in fi. It's Laplace transform 

(2.4) H g , e (wi,...,w £ ) = fi 5 (M)e-^^ 1+1 H™+w(« < +i)) 

is the function we consider in this paper. We note that the automorphism group Aut(p) 
acts trivially on the function e _ ^ 1 ( Wl+1 ) H w j 1 i c j 1 explains its appearance in (|2.3[) . 
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Since the coordinate change x = e~ ( - w+1 ^ identifies 

(2.5) -Q——^H g ^(w(x 1 ),...,w(x e )) = H gie (x 1 , . . . ,x e ), 

the Laplace transform (12.41) is a primitive of the generating function (|1.4j) . 

Before performing the exact calculation of the holomorphic function TC g .e(wi, . . . , we), let 
us make a quick estimate here. From Stirling's formula 

uk+n i 

e fe — = k 2, 

k\ V2^ 

it is obvious that 7i g ^(wi, . . . , we) is holomorphic on Re(wi) > for all i = 1, . . . ,£. Because 
of the half-integer powers of /Vs, the Laplace transform Ti. g ^(wi, . . . , we) is expected to be a 
meromorphic function on a double-sheeted covering of the u>j-planes. Such a double covering 
is provided by the Lambert curve C of (jl.3p . So we define 

OO j 

(2.6) t = l + Eji e - k(W+1) > 

k=l 

which gives a global coordinate of C. The summation converges for Re(w) > 0, and the 
Lambert curve is expressed in terms of w and t coordinates as 

(2.7) e -(-+i) = ^i_^ e -(l-|). 

The w-projection tt : C — > C is locally a double-sheeted covering at t = oo. The inverse 
function of (|2.6p is given by 

(2.8) w = w(t) = -l-l 0g ( 1 -l) = f : li., 

v 7 m=2 

which is holomorphic on Re(t) > 1. When considered as a functional equation, (|2.8|) has 
exactly two solutions: f and 

, . , , 2 4 r> 8 O 8 A 

(2.9) s(t) = -tH 1 i~ 2 H t" 3 H i" 4 H . 

v ; w 3 135 405 567 

This is the deck-transformation of the projection tt : C — > C near t = oo and satisfies the 
involution equation s(s(i)) = t. It is analytic on C \ [0, 1] and has logarithmic singularities 
at and 1. Although w(t) = w[s(t)), s(t) is not given by the Laplace transform (12. 6j) . 
Since the Laplace transform 

00 yk+n 

(2.10) | n (t) = ^ ^rpe-^+l) 

fe=l 

also naturally lives on C, it is a meromorphic function in t rather than in w. Actually it is 
a polynomial of degree 2n + 1 for n > because of the recursion formula 

(2.11) L+l{t) =t 2 (t-l)- Ut) foralln>0, 

at 

which follows from ([2TTU]) . and d^SJ that implies 
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We note that the differential operator of (|2.1ip is discovered in [T7]. For future convenience, 
we define 

(2-13) i_ 1 ( t ) = t -^=y, 

which is indeed the y coordinate of the original Lambert curve (jl.3p . We now see that the 
Laplace transform 

(2.i4) n 9/ (h, ...,k) = n g , t {(w(ti), . . .,w(k)) = E ^^)e~^(- i+1 ) + - + «^ +1 )) 

n\-\ Vn t <3g-?,+l i=l 

is a symmetric polynomial in the t- variables when 2g — 2 + £ > 0. 

It has been noted in [Tj I8| I§1 \TT\ I12j that the Airy curve w = ^v 2 is a universal object of 
the topological recursion for the case of a genus spectral curve with only one critical point. 
Analysis of the Airy curve provides a good control of the topological recursion formula for 
such cases. The Airy curve expression is also valid around any non-degenerate critical point 
of a general spectral curve. To switch to the local Airy curve coordinate, we define 

function in t that solves 

/ \ 1 2 1 , A 1\ 1 , / 1 \ A 1 1 

(2.16) - v 2 = w = log 1 - - = - log 1 - — = > . 

V ^ 2 t *\ t) s(t) 6 V *(t)J ^ 2 ™t™ 

Note that we are making a choice of the branch of the square root of w that is consistent 
with (12. 6p . The involution (12. 9ft becomes simply 

(2.17) v{t) = -v{s(t)). 

The new coordinate v plays a key role when we later reduce the Laplace transform of the 
cut-and-join equation (|1.2p to the Bouchard-Marino topological recursion (jl.5p . 

3. The cut-and-join equation and its Laplace transform 

The cut-and-join equation for Hurwitz numbers was discovered in |151I37|. It has become 
an effective tool for studying algebraic geometry of Hurwitz numbers and many related 
subjects [HI [T71 [211 [2H ESI [39]. In this section we calculate the Laplace transform of the 
cut-and-join equation following [32J, and prove Theorem ll.il 

The simplest way of presenting the cut-and-join equation is to use a different primitive 
of the same generating function of Hurwitz numbers (|1.4|) . Let 

(3.1) H(s,p) = EE H 9< £ ( S ' P )' H 9,H S >P)= E h 9,nPt* — > 

g>o e>i n-x(n)=i 

where = PuxPu-z " " 'Puei anc l r = 2^ — 2 + £+ |^| is the number of simple ramification 
points on P. The summation is over all partitions of length £. Here pk is the power-sum 
symmetric function 

(3.2) Pfe = E^' 

i>l 
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which is related to the monomial symmetric functions by 



TeS e i=l 

Therefore, we have 

— -^-H 9 ,(l,p) = ^(x 1 ,...,x,)= J; (2 g - 1 2 + y + H)! ^^n<5 1 ' 

which is the generating function of (|1.4p . Because of the identification (|2.5p . the primitives 
H Sj ^(l,p) and "Hg^itx, . . . ,t() of (|2.14j) are essentially the same function, different only by 
a constant. 

Remark 3.1. Although we do not use the fact here, we note that H(s, p) is a one-parameter 
family of r-functions of the KP equations with ^p^ as the KP time variables \21\ [M] . The 
parameter s is the deformation parameter. 

Let z £ P 1 be a point at which the covering / : X — > P 1 is simply ramified. Locally we 
can name sheets, so we assume sheets a and b are ramified over z. When we merge z to oo, 
one of the two things happen: 

(1) The cut case. If both sheets are ramified at the same point Xi of the inverse image 
f~ 1 (oo) = {x\, . . . ,X(}, then the resulting ramification after merging z to oo has a 
profile 

(jii, . . . . . . ,n e ,a,m - a) = (/j,(i),a,m - a) 

for 1 < a < Hi. 

(2) Otherwise we are in the join case. If sheets a and b are ramified at two distinct 
points, say Xi and Xj above oo, then the result of merging creates a new profile 

(fix, ...,£&,... ,ftj, ...,fii,m + fij) = + fij). 

Here the ~ sign means removing the entry. The above consideration tells us what happens 
to the generating function of the Hurwitz numbers when we differentiate it by s, because 
it decreases the degree in s, or the number of simple ramification points, by 1. Since the 
cut case may cause a disconnected covering, let us use the generating function of Hurwitz 
numbers allowing disconnected curves to cover P . Then the cut-and-join equation takes 
the following simple form: 



d_ 

ds 2 



a,/3>l 

It immediately implies 



^E ((« + ^^ + *-^) 



e 



H(a,p) _ Q 



(3 - 3) *r = 2 E r + + "fo** + a ^dp~ a " dp-J ' 

which is the cut-and-join equation for the generating function H(s,p) of the number of 
connected Hurwitz coverings. 

Let us now apply the ELSV formula (|1.1|) to (|3.1|) . We obtain 



1 _ _ ..HT". 

(3.4) H 9/(s , P )=i y: (^A:(i)) 9i / 9 -^nEV s ^ 



LAPLACE TRANSFORM OF THE CUT-AND-JOIN EQUATION 



i £ HMW T = E ' 



where H g (fj,) is introduced in (|2.3|) . Now for every choice of r > 1 and a partition //, the 
coefficient of p^s 7 * -1 of the cut-and-join equation (|3.3p gives 

Theorem 3.2 ([32]). The functions H g ((j,) of (|2.3p satisfy a recursion equation 
(3.5) r{g,n)H g (ji) = ^2(fH + fj, j )H g (fx(i,j),iJ, i + ij, j ) 



+ 9 E E ^ 



i=l a+/3=m 

Remark 3.3. Note that 



\ 



H g -i{n(i),a,0) + E H giiyi^)Hg 2 (v 2 ,P) 



91+92=9 



£{fx(i,j)) =£-2 

£(u 1 )+£{u 2 )=£( f x(i))=l-l. 

Thus the complexity 2g — 2 + £ is one less for the coverings appearing in the RHS of 
(I3.5|) . which is the effect of applied to H(s,p), except for the unstable geometry 

corresponding to = and = in the join terms. If we move the (0, l)-terms to the 
LHS, then the cut-and-join equation (13. 5p becomes a topological recursion formula. 

Let us first calculate the Laplace transform of the cut-and-join equation for the £ = 1 
case to see what is involved. We then move on to the more general case later, following 

Proposition 3.4. The Laplace transform of the cut-and-join equation for the 
gives the following equation: 

(3.6) Yl ^ g {l)) 9tl [(2g-l)Ut)+L+i(t){l-Li(t)) 

n<3g-2 

stable 



1 case 



2 ^ 

a+b<3g-4 



(r^A^l)}^ + E <^A^(l)) Sljl <T- b A^(l)) 52;1 

91+92=9 



£a+l(i)6+l(*)- 



Proof. The cut-and-join equation for £ = 1 is a simple equation 
(3., 



.7) (25-1 + ^(^ = 1 £ (a, /?)+ E H gi (a)H g M) ■ 

a+/3=n \ 91+92=9 ) 



The Laplace transform of the LHS of (|3.7|) is 

E <7nA, v (l)> 9>1 [(2g - l)Ut) + L+i(t) ■ 

n<3g-2 

When summing over /j, to compute the Laplace transform of the RHS, we switch to sum 
over a and (5 independently. The factor | cancels the double count on the diagonal. Thus 
the Laplace transform of the stable geometries of the RHS is 

stable 



2 ^ 

a+b<3g-i 



(T 8 7 k AjL 1 (l)) Jh . 1)2 + E ^Al(l)) gi>1 (nA v g2 (l)) 92il 

91+92=9 



e a+ i(i)6+i(i). 



12 B. EYNARD, M. MULASE, AND B. SAFNUK 

The unstable terms contained in the second summand of the RHS of (|3.T|) are the g = 
terms Ho(a)H g ((3) + H g (a)Ho(f3). We calculate the Laplace transform of these unstable 
terms using (|2.ip . Since 

a a ~ 2 

Ho(a) = — r , 
a! 

the result is 

J2(r a ^ g (l)) gil U(t)i a+1 (t). 

a 

This completes the proof. □ 

Remark 3.5. We note that (|3.6p is a polynomial equation of degree 2n + 2. Since C~i(t) = 
1 — t, the leading term of £ n +i(i) is canceled in the formula. 

To calculate the Laplace transform of the general case (|3.5|) . we need to deal with both of 
the unstable geometries (g,£) = (0,1) and (0,2). These are the exceptions for the general 
formula ([2TT3]) . Recall the (0, 1) case ([23]) . The formula 

00 l.fc-2 i 

(3.8) Wo,i(*) = E "Tr^^ = -^72 + c = 

fc=i 

where the constant c is given by 

00 j.fc-2 

fc=i 

is used in ([3.6H . The (g,^) = (0,2) terms require a more careful computation. We shall see 
that these are the terms that exactly correspond to the terms involving the Bergman kernel 
in the Bouchard-Marino recursion. 

Proposition 3.6. We have the following Laplace transformation formula: 

.Ml ,,M2 



(3.9) H„2(*i,*2)= F — ^ • • ^ 7 e-w(^ +1 ) e -«'(^+ 1 ) 
^ + Ml' ^2! 



^ g-l(tl)-g-l(t 2 ) \ i - 

= log - - 4-l(*2>- 

Proof. Since x = e"^ 4 " 1 ^, (j3.9p is equivalent to 

1 ,,M2 / 00 J.fc-1 fc k\ 

(3.10) £ — L_ •^ T -^xf=log ' 

(^1,^2)^(0,0) 

where \xi\ < e _1 , |x2| < e" 1 , and < \x\ — X2\ < e~ 1 so that the formula is an equation of 
holomorphic functions in x\ and xi- Define 

( M i,M2)^(0,0) 

Then 



,o)^E^.-io g 

/Ltl>l ^ \fc=l 
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l-i(t) - log (^^) = 1 - \ ~ log (l - jj + logx 

1 / V 
= 1 log 1 I — w 



t V t 

due to (|2.8p . Here t is restricted on the domain Re(t) > 1. Since 




x 



x\ - x 2 



- 1)^ ^g (1-1 (tl) - 1-1 (*2)) " XI— l0g(xi - X 2 ) 



ti(*i-l)^-log -- + - 



<9ii V *i *2 / ^1 — ^2 



ti — t 2 £1 — ^2 



we have 



x l a !" — 




X\ - x 2 



t\ —t 2 Xl — X2 

= ht 2 - 1 = ia(h)io(t 2 ) + fo(ii) + Uh)- 



On the other hand, we also have 

dx\ dx 2 ^ M1+M2 Mi' M2! 



(Ml,M2)^(0,0) 



&(*l)6(*2)+6(*l)+&(*2)- 



Therefore, 

(3.1!) (, 1 ^_ +I2 ^_) 0(X1 , X2) = O . 

Note that </>(xi,X2) is a holomorphic function in x\ and cc2- Therefore, it has a series 
expansion in homogeneous polynomials around (0,0). Since a homogeneous polynomial in 
xi and x 2 of degree n is an eigenvector of the differential operator X\-Mr + ^2gf- belonging 
to the eigenvalue n, the only holomorphic solution to the Euler differential equation (|3.1ip 
is a constant. But since 4>(xi,0) = 0, we conclude that (f>{x\,X2) = 0. This completes the 
proof of (|3.10p . and hence the proposition. □ 



The following polynomial recursion formula was established in [32]. Since each of the 
polynomials TO(£l)'s in (|3. 12|) satisfies the stability condition 2g — 2 + 1 > 0, it is equivalent 
to (jl.2p after expanding the generating functions using (|2,14p . 

Theorem 3.7 (|32j). The Laplace transform of the cut- and- join equation (|3.5p produces 
the following polynomial equation on the Lambert curve: 
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(3.12) (^g - 2 + i + i 1 ~ Li{U))tUU ~ l)^: J Hgjfa, ...,t e ) 

ti(ti ~ l) 2 ^7^9/-i (*!>•••)*/)•••! tl) — tj(tj — l) 2 ^-H g /-i (h, . . . ,ti, . . . ,te) 



i<j 



a 



y^^fiu - 1)^7^5,^-1 (*i> ■ ■ ■ ,tj, . . . ,t e ) 



i=l 



i - r d 2 ^ 

+ o Z, «l(«l - l)«i(«2 - l) o o ftg-M+l "2, *L\{<}) 



du\du2 



Ul=U2=U 



stable 







i=l 31+32=9 
JUi^=L\{i} 



St- 



In t/ie last sum eaca term is restricted to satisfy the stability conditions 2g± — 1 + | J| > 
and 2a 2 - 1 + |if | > 0. 

Remark 3.8. The polynomial equation (|3.12p is equivalent to the original cut-and-join 
equation (|3.5[) . Note that the topological recursion structure of (|3.12p is exactly the same 
as (jl.5p . Although (|3.12p contains more terms, all functions involved are polynomials that 
are easy to calculate from (|2.1ip . whereas (|1.5|) requires computation of the involution s(t) 
of (|2.9p and infinite series expansions. 

Remark 3.9. It is an easy task to deduce the Witten-Kontsevich theorem, i.e., the Virasoro 
constraint condition for the ^-class intersection numbers [38^ I22j. from (|3.12p . Let us use 
the normalized notation a n = (2n + l)!!r n for the "0-class intersections. Then the formula 
according to Dijkgraaf, Verlinde and Verlinde (see for example, [1]) is 



(3.13) (a n ar, 



nL/g,e+l 



n L\{i}lg,l 



a+b=n—2 



\ Y l \ a ci<y b a nL ) g _ 1/+2 + Y{2ni + l){(T n+ni _n 

ieL 

- stable 



31+92=3 a+b=n—2 
I\JJ=L 



Eqn. (|3.13p is exactly the relation of the homogeneous top degree terms of (|3.12p . after 
canceling the highest degree terms coming from £ ni+ i(tj) in the LHS 



4. The Bouchard-Marino recursion formula for Hurwitz numbers 

In this section we present the precise statement of the Bouchard-Marino conjecture on 
Hurwitz numbers. 

Recall the function we introduced in (12.61): 



(4.1) 



k=l 



k k 

¥' 



This is closely related to the Lambert W -function 



(4.2) 



W(x) 



k=l 
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By abuse of terminology, we also call the function t(x) of (|4,ip the Lambert function. The 
power series (I4.ip has the radius of convergence 1/e, and its inverse function is given by 

(4.3) x = x (t) = ^ Tl-ijei 

Motivated by the Lambert VF-function, a plane analytic curve 

(4.4) C = {(x, t)\x = x(t)} cCxC 

is introduced in [3], which is exactly the Lambert curve (|1.3p . We denote by it : C — > C the 
x-projection. Bouchard-Marino [3] then defines a tower of polynomial differentials on the 
Lambert curve C by 

(4.5) Ut) = ~ 



t 2 (t-l)i n ^{t) 



with the initial condition 

(4.6) Ut) = dt. 

It is obvious from (|4,5p and (|4.6p that for n > 0, is a polynomial 1-form of degree 2n 

with a general expression 



Ut) = t r ' 



(2n + 1)!! t n - ( ( 2re + 3 ) !! - (2n + 1)!! ) f 1 " 1 + • • • + (-l) n (n + 1)! 



dt. 



All the coefficients of £ n (t) have a combinatorial meaning called the second order reciprocal 
Stirling numbers. As we will note below, the leading coefficient is responsible for the Witten- 
Kontsevich theorem on the cotangent class intersections, and the lowest coefficient is related 
to the A 5 -formula. For a convenience, we also use = t~ 2 dt and £-2(2) = t~ 3 dt. 

Remark 4.1. The polynomial £ n (t) of (|2.10p is a primitive of £ n (i): 

(4.7) d£„(i) = Ut)- 

Definition 4.2. Let us call the symmetric polynomial differential form 

e 

d® i n g , e (ti,...,u)= Yl fou---vW)> C8>e«*(*<) 

rtiH h»i£<3g— 3+f i=l 

on C £ the Hurwitz differential of type (g,£). 

Remark 4.3. Our £ n (t) is exactly the same as the £ n (y)-differential of [3]. However, this 
mere coordinate change happens to be essential. Indeed, the fact that our expression is a 
polynomial in t- variables allows us to calculate the residues in the Bouchard-Marino formula 
in Section [5j 



Remark 4.4. The degree of d® 7i g /(ti, . . . , ti) is 2(3<? — 3 + £), and the homogeneous top 
degree terms give a generating function of the ^-class intersection numbers 

e i 

< r »i • • • T ^> II( 2ni + ^ ® dt i- 

niH \-n e =3g-3+£ i=l 1=1 

The homogeneous lowest degree terms of d® n. g ,e(ti, . . . ,ti) are 

i i 
( _ 1)3ff -3+^ £ (r ni ---r n ,A 9 ) Hini + iy.tr^dU. 

niH hn^=2g-3+£ i=l i=l 
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The combinatorial coefficients of the A 9 -formula [13} Q3] can be directly deduced from the 
topological recursion formula (|1.2p |32j . 

Remark 4.5. The unstable Hurwitz differentials follow from (|2.ip and ()3.9p . They are 
(4.8) dH QA (t) = ^ dt = t-2(t); 

{ti-t 2 y {xi-x 2 y 

Remark 4.6. The simplest stable Hurwitz differentials are given by 

d® z HQ,z{ti,t 2 ,h) = dh ® dt 2 ® dt 3 ; 
(4-10) ^ i i 

dHi,!(t) = - ( - £o(t) + £(*)) = 24 (* - l)(3t + 1) dt. 

The amazing insight of Bouchard and Marino [3] is that the Hurwitz differentials of 
Definition 14.21 should satisfy the topological recursion relation of Eynard and Orantin [11] 
based on the analytic curve C of (14. 4p as the spectral curve. Since the topological recursion 
utilizes the critical behavior of the x-projection tt : C — > C*, let us examine the local 
structure of C around its critical points. Let z = — j be a coordinate of C centered at 
t = oo. The Lambert curve is then given by 

x = -(1 + z)e~ z . 
e 

We see that the x-projection tt : C — > C* has a unique critical point qo at z = 0. Locally 
around qo the curve C is a double cover of C branched at go- For a point q £ C near qo, let 
us denote by q the Galois conjugate point on C that has the same x-coordinate. Let S(z) 
be the local deck-transformation of the covering tt : C — > C*. Its defining equation 

(4.11) 5(2) - log (1 + 5(z)) = z - log(l + z) = y y -^z m 

' L — ' TO 

m=2 

has a unique analytic solution other than z itself, which has a branch cut along (— oo, — 1]. 
We note that S is an involution 5(5(z)) = z, and has a Taylor expansion 

2 2 4 „ 44 4 104 - 40 B 7648 , 2848 s ^. ,. 

S(z) = -z + -z 2 - -z 3 H z 4 z 5 H z 6 z 7 H z 8 + 0(z 9 ) 

w 3 9 135 405 189 42525 18225 v ' 

for \z\ < 1. In terms of the i-coordinate, the involution corresponds to s(t) of f) 2 . 9 1) : 

/*(«) = -i=* 

The equation (|4.1ip defining 5(z) translates into a relation 

(4.12) = - . = -du, = -t,cfo = tt* f^V 

v ; * 2 (i-i) s (t) 2 ( s (t)-i) w 

Using the global coordinate t of the Lambert curve C, the Bergman kernel is defined by 
4.13 S(ti,t 2 ) = tt , si =dtAlog ti-i2 . 

We have already encountered it in (|4.9p in the expression of Tto, 2 (ti,t 2 ). Following [TT] . 
define a 1-form on C by 
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dE ^ ^) = \[ q B(-M = \ - dt 2 

= Ui-i( s (h) - 1 2 ) -i-ifa-pAdh, 

where the integral is taken with respect to the first variable of B(ti,t 2 ) along any path from 
q to q. The natural holomorphic symplectic form on C* x C* is given by 

O = d log y A d log x = d log ^1 — ^\ Ad log x. 

Again following O [TT] , let us introduce another 1-form on the curve C by 

"<«•* - jf n < ' ■«<«» - - 4) wih) = («-( s < ( ») " «-'< ( >)^T)' 

The kernel operator is defined as the quotient 

dE(q,q,t 2 ) 1 ti s(h) gfti - 1) 

K{ti,t 2 ) = 7 — =r — = - • r— — — ® dt 2 , 

Lo{q,q) 2 ti-t 2 s{t\) — t 2 dh 

which is a linear algebraic operator acting on symmetric differential forms on C by replacing 
dti with dt 2 . We note that 

(4.14) K(t 1 ,t 2 ) = K(s(t 1 ),t 2 ) , 

which follows from (|4.12p . In the z-coordinate, the kernel has the expression 

(4.15) K = ---^^--. r-f- — -.dt®l- 

K J 2 z (l + zt)(l + S(z)t) dz 

1/1 1, x 1/ x 9 

= - J ( - + 1 + 3 (3* " 2)tz + -(3t - 2)tz 2 

+ — (135t 3 - 180i 2 + m + 16)tz 3 + • • • ] dt ® — . 
135 / dz 

Definition 4.7. The topological recursion formula is an inductive mechanism of defining a 
symmetric £-form 

W g/ (t L ) = W 9te (t 1 ,...,t e ) 
on C e for any given g and ^ subject to 2g — 2 + £ > by 



(4.16) W g/+l {t ,t L ) = -^-j 



K(t,t )\W g _ l/+2 (t,s(t),t L ) 

i=l 

stable terms x 

+ E %i.[i[+i(t.*/)®w r fl2 ,|j[+i(a(*).*j)) 

31+92=9, /LI 

Here ij = (ii)ie/ for a subset I C L = {1,2,... ,£}, and the last sum is taken over all 
partitions of g and disjoint decompositions I\JJ = L subject to the stability condition 



7oo 
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2<7i — 1 + \I\ > and 2g 2 — 1 + | J\ > 0. The integration is taken with respect to dt on the 
contour 700, which is a positively oriented loop in the complex t-plane of large radius so 
that \t\ > max(|i |, |s(io)l) f° r * ^ 7oo- 

Now we can state the Bouchard-Marino conjecture, which we prove in Section [7J 

Conjecture 4.8 (Bouchard-Marino Conjecture [3]). For every g and I subject to the sta- 
bility condition 2g — 2 + I > 0, the topological recursion formula (|4.16p with the initial 
condition 

(4.17) 

gives the Hurwitz differential 



W ,3(ti,t 2 ,t 3 ) = dh ® dt 2 <g> dt 3 
Wi,i(h) = ifa- 1) (3ti + l) dh 



W g , e (t 1 ,...,t e ) = d® e H g ,e(t 1 ,...,t e ). 
Remark 4.9. In the literature OUT], the topological recursion is written as 

dE(q,q, t ) 



(4.18) W g/+ i(t ,t L ) =Res q=q 



u>(q,q) 



W g „ 1/+2 (t(q),t{q),t L ) 



+ J2 W gumi (t{q)M)®W g2 ^ +1 (t(q),tj) 

gi+g2=g, l\JJ=L , 
including all possible terms in the second line, with the initial condition 

|Vo,i(ii) = 
lWo )2 (ii,£ 2 ) = B(tx,t 2 ). 



(4.19) 



If we single out the unstable terms from (|4.18p . then we obtain (|4.16p . Although the initial 
values of W g ^ given in (|4.19|) are different from (|4.8p and (|4.9p . the advantage of ()4.18p is 
to be able to include (|4,17p as a consequence of the recursion. 

Remark 4.10. It is established in [11] that a solution of the topological recursion is a 
symmetric differential form in general. In our case, the RHS of the recursion formula (|4.16p 
does not appear to be symmetric in to,t\, . . . ,te. We note that our proof of the formula 
establishes this symmetry because the Hurwitz differential is a symmetric polynomial. This 
situation is again strikingly similar to the Mizrakhani recursion [29, 30J, where the symmetry 
appears not as a consequence of the recursion, but rather as the geometric nature of the 
quantity the recursion calculates, namely, the Weil-Petersson volume of the moduli space 
of bordered hyperbolic surfaces. 

5. Residue calculation 

In this section we calculate the residues appearing in the recursion formula (|4.16p . It 
turns out to be equivalent to the direct image operation with respect to the projection 
vr:C7^C. 

Recall that the kernel K(t,to) is a rational expression in terms of t,s(t) and to. The 
function s(t) is an involution s(s(t)) = t defined outside of the slit [0,1] of the complex 
t-plane, with logarithmic singularities at and 1. Our idea of computing the residue is to 
decompose the integration over the loop 700 into the sum of integrations over 700 — 7[o,i] 
and 7[o,i]> where 7[o,i] is a positively oriented thin loop containing the interval [0, 1]. 
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Figure 5.1. The contours of integration. 700 is the circle of a large radius, and 
7[o,i] is the thin loop surrounding the closed interval [0, 1]. 

Definition 5.1. For a Laurent series ^nez a ^ n ' we denote 



n>0 



Theorem 5.2. In terms of the primitives £ n (t), we have 



(5.1) R afi (t) 



2-ni 



K(t',t)Ut'Ms(t')) 



7oo 



ts{t) 
t - s(t) 



£a(t)6+l (s(t)) + L+l («(*))&(*) 



Similarly, we have 

(5.2) R n (t,t i ) = - — I K(t', t) (^ n (t')B(s(t'), ti) + £ n (s(t'))B(t', U) 



ts(t) 
t - s(t) 



£n+l 



Proof. In terms of the original z-coordinate of [3] , the residue R a ^ (t) is simply the coefficient 
of z^ 1 in K{t' ,t)^ a (t')^b(s(t')) , after expanding it in the Laurent series in z. Since £«(£') is 
a polynomial in t' = —-, the contribution to the z _1 term in the expression is a polynomial 
in t because of the z-expansion formula (|4.15p for the kernel K. Thus we know that R a ,b{t) 
is a polynomial in t. 

Let us write £, n (t) = f n (t)dt, and let 7[o,i] be a positively oriented loop containing the slit 
[0, 1], as in Figure 1531 On this compact set we have a bound 



ts(t) „ 



t 2 (t-l)s'(t)f a (t)f b (s(t)) 



< M, 



t - s{t) 

since the function is holomorphic outside [0, 1]. Choose |t| >> 1. Then we have 



K(t',t)Ut'Ms(t')) 

Z7Tl J 1[0,1] 
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1 

2vri 



< 



M 
2^ 



7[0,l] 



7[o,i] 

1 

2 



1 /_1 t'sjt') 

2 \ f - t ~ 8 (f) - t ) s(t') - t 



7 t' 2 (t'-l)s'(t')f a (t')f b (s(t'))dt' 



(g> dt 



1 



1 



t' - t s(t') - t 



,/ , M , 

dt (8) dt ~ — — 
4?r t 



Therefore, 
1 

~2tH 



I K(t',t)Ut'Ms(t')) = -^~i K(t',t)Ut'Ms(t'))+0(t- 1 

Too Z7ri ^7oo-7[0,ll 



1 

"2vri 



•/ Too— Trn n 



'Too— 7[o,i] 

Noticing the relation (|4.12p and the fact that s(t) is an involution, we obtain 
1 

~ 2tH 



K(t',t)Ut'Ms(t')) 



Too — f[0,l] 



1 

2tt7 

1 

27d 



1 



t's(t') 



-( — - 



7 t' 2 (t'-l)s>(t>)f a (t')f b (s(t'))dt'®dt 



Too-T[0,l] 



+ 



1 

2^i 



1 1 t's(t') 
x] 2 ' t'-t ' s(t') - t' 

1 1 t's(t') 



1 tg(t) „ 2 

2 ' t - s(t) 



(7<»)-(7 [ o.i]) 2 'W"* S (*')-* 
•t 2 (t-iy(t)/ a (t)/ fc ( S (t))dt 



t ,2 (t' - l)s'(t')f a (t')f b (s(t'))dt' ® dt 



r 1?tf-l)fatf)fb(8tf))d8(lf)®dt 



+ * • • ^(t) 2 (sW - l)/a(s(t))/ 6 (i)dt 



2 t - s(t) 
* s O t 2 (t-l) s '(t) 



fa(t)fb(s(t)) + f a (s(t))f b (t) 



t - s(t) 
1 ts(t) 



dt 



2 t - s(t) 

1 ts(t) 

2 ' t - s(t) 

1 ts(t) 

2 ' t - s(t) 



e (*)Cfr4-i («(«)) («(«))&(*) 



r/t 



| a+1 (t)4+l («(*)) + |a+l (s(t))6+l(*) 



dt 



Here we used (|2.1ip and f)4. 12j) at the last step. The proof of the second residue formula is 
exactly the same. □ 

Remark 5.3. The equation for the kernel (|4,14p implies 

Ra,b(t) = R bAt) = ~ [Ra,b{a(t))] + . 

Let us define polynomials P ab (t) and P n (t,ti) by 



(5.3) P a ,b(t)dt 



dt 



ts(t) 
t - s(t) t 2 (t - 1) 



L+i(t)ib+i{s(t))+L+i{s(t))i b+1 (t) 
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(5.4) P n (t,U)dt®dti = dt 



ts{t) ( £ n+1 (t)ds(t) i n+ i(s(t))dt 



t - s(t) \ s(t) -U t-U 



Obviously deg P a ,b(t) = 2(a + b + 2). To calculate P n (t,U), we use the Laurent series 
expansion 



oo /, \ k 



t-U t 

k=0 



and take the polynomial part in t. We note that it is automatically a polynomial in U as 
well. We thus see that deg P n (t, U) = 2n + 2 in each variable. 

Theorem 5.4. The topological recursion formula (|4.16j) is equivalent to the following equa- 
tion of symmetric differential forms in I + 1 variables with polynomial coefficients: 



Y ( T nTn L A V g (l)) g/+1 dUt)®d£ nL (t L ) 



n,n L 



= Y Y (VnAW^l^sA^^)*®*®^^^}) 

+ ( Y (^b^L^-li 1 )) g-l,£+2 
\ a,b,riL 

stable \ 

+ Y Y (^^.(^./.j/i+iK^i^toM+i \p a ,b{t)dt ® ^(ti). 

gi+92=g / 
.Here L = {1, 2 . . . ,£} is an index set, and for a subset I C L, we denote 

ti = (ti) ie i, ni = {ni\iel}, T ni =Y\_T nt , d£ ni (ti) = (^-^-( ni (U)dU- 

The last summation in the formula is taken over all partitions of g and decompositions of 
L into disjoint subsets I\\J = L subject to the stability condition 2g\ — 1 + \I\ > and 
2<? 2 -l + |J| >0. 

Remark 5.5. An immediate observation we can make from (|1.5p is the simple form of the 
formula for the case with one marked point: 

(5.6) ^ g iX)) 9>1 j t Ut) 

n<3g-2 

stable 

= Y ((^ A ^l( 1 )) 9 -l,2+ Y (^ 1 (l)) fllil (7 i AV(l)) flail )p a(6 (t). 
a+fe<3g-4 9i+92=g 

6. Analysis of Laplace transforms on the Lambert curve 

As a preparation for Section [7] where we give a proof of (jl.5p , in this section we present 
analysis tools that provide the relation among the Laplace transforms on the Lambert curve 
(|2.7p . The mystery of the work of Bouchard-Marino [3] lies in their £ n (y)-forms that play 
an effective role in devising the topological recursion for the Hurwitz numbers. We have 
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already identified these differential forms as polynomial forms d£ n (t), where £ n (i)'s are the 
Lambert W-function and its derivatives. 
Recall Stirling's formula 

(6.1) logr(z) = ^log27r+r«-^logz-* 

A B 2r 2r+1 _ J_ [°° B 2m (x-[x]) 
Z_, 2r(2r - 1) 2m 7 (2 + x) 2m 

where m is an arbitrary cut-off parameter, B r (s) is the Bernoulli polynomial defined by 

ze zx A , , 2 r 

— =>B r (x) — , 

e z — 1 r! 

B r = B r (0) is the Bernoulli number, and [x] is the largest integer not exceeding x £ R. For 
N > 0, we have 



\/2^ V 4t2r(2r-l) i \2m J (N + x) 2m J 



■f 2r(2r 



Let us define the coefficients for fc > by 

00 / 00 „ 

( 6 .2) 5>*r* = «P -yj-^-iv- 



fc=0 



2r(2r - 1) 



12 288 51840 2488320 

Then for a large iV we have an asymptotic expansion 

N N+n -I 00 

(6.3) ^^—.-L^-lg^. 
Definition 6.1. Let us introduce an infinite sequence of Laurent series 

k=0 



for every n G Z, where s^'s are the coefficients defined in (|6,2p . 

The following lemma relates the polynomial forms £ n (t) = d£ n (t), the functions r) n {v), 
and the Laplace transform. 

Proposition 6.2. For n > 0, we have 

roc 

(6.5) / e~ s — —e~ sw ds = rj n (v) + const + 0(w) 

Jo r (s + l) 

with respect to the choice of the branch of y/w specified by v = —y2w as in (|2.8[) and 
(|2.15p . where 0(w) denotes a holomorphic function inw = ^v 2 defined around w = which 
vanishes at w = 0. The substitution of (|2.15p in %(i>) yields 

(6.6) n» = i U» - %(-^)J = \ Un(t) - L{s(t)) 



(6.7) r ] ^(v) = UU(t)-U{s(t) 
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where s(t) is the involution of (|2.9p . This formula is valid for n > — 1, and in particular, 
we have a relation between the kernel and rj-i(v): 

More precisely, for n > —1, we have 

'Vn(v) =in(t)+F n (w) 
K Vn(-v) =L{s(t)) +F n {w) ' 

where F n (w) is a holomorphic function in w. 

Proof. From definition (|6.4p . it is obvious that the series rj n (v) satisfies the recursion relation 

(6.9) Vn+i(v) = -- -r-rintv) 

v dv 



(6. 



for all n E Z. The integral (16. 5p also satisfies the same recursion for n > 0. So choose an 
n > 0. We have an estimate 



Q s+n I 

^ s "-^5> s - fc + o(^ 



r(. + i) v^f A ._ u 



2 



that is valid for s > 1. Since the integral 



1 „s+n 

e~ s -^ -e~ sw ds 

r( s + i 





is an entire function in w, we have 

oo „s+n pi „s+n poo „s+n 

e~ s -^ -e- sw ds= / e~ s -^ -e- sw ds+ / e~ s -j- -e~ sw ds 

r(s + i) Jo r( s + i) A r( 8 + i) 

= / — s"^Vs feS -Me- M (ls+ / 0(s - 2)e~ sw cfa + const + 0(w) 

JO Vv27T ^ / Jl 

1 » (2(n-fc)-l)!! 
= - 2^ Sfc + const + uO(w) + O(w). 

This formula is valid for all n > 0. Starting it from a large n » and using the recursion 
(|6.9p backwards, we conclude that the vO(w) terms in the above formula are indeed the 
positive power terms of rj n (v). The principal part of r] n {v) does not depend on the addition 
of positive power terms inw= \v 2 , since — transforms a positive even power of v to a 
non-negative even power and does not create any negative powers. This proves (|6.5p . 

Next let us estimate the holomorphic error term 0{w) in (|6.5p . When n < —1, the 
Laplace transform (|6.5p does not converge. However, the truncated integral 

„s+ra 

~ s -^ -e~ sw ds 

T(s + 1) 

always converges and defines a holomorphic function in v = —y/2w, which still satisfies the 
recursion relation (16. 9p . Again by the inverse induction, we have 

poo „s+n 

(6.10) / e- s ——-e-™ds = Vn (v) + 0(w) 

Ji r(s + i) 

for every n < 0. Now by the Euler summation formula, for n < — 1 and Re{w) > 0, we have 
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roc 



„s+n 00 uk+n 

= -\<>~ {w+1) + /°° (' ~ I s ' ~ ^) S ( e "r^Ti) <i ~™') 

Note that the RHS of (|6.1ip is holomorphic in w around w = 0. Prom (|2.10p . (|6.10p and 
(|6.1ip . we establish a comparison formula 

(6.12) r j - 1 (v)-€- 1 (t)=F-. 1 (w), 

where F-i(w) is a holomorphic function in w defined near w = 0, and we identify the 
coordinates t,v and w by the relations (|2.16p and (|2.15p . Note that the relation (|2.16p is 
invariant under the involution 



(6.13) 

Therefore, we also have 
Thus we obtain 



v 1 — > —v 
t I — ► s(t) 



V-l(v) = Ui-i(t)-£-i{s(t)) 



which proves (|6.7p . Since — = t 2 {t — the recursion relations (|6.9p and (|2.1ip for 
£, n (t) are exactly the same. We note that from (|4.12p we have 

' t*(t-l)± = S (tf(s(t)-l)- 



v dv dt ds(t) 

Therefore, the difference £n+i(t) — ^ n +i(s(t)) satisfy the same recursion 

(6.14) £ n+1 (t) - i n+1 (s(t)) =t\t-l)j(Ut)-£ n (s(t))\ 

The recursions (I6.9jl and (|6.14p . together with the initial condition (|6.7p . establish (|6.6p . 
Application of the differential operator 

Id d 9/ . d 

r = — r- =* *-!)-r 

v at> ait; at 

(n + l)-times to (|6.12p yields 

%(«) - Lit) = F n (w), 

where F n (w) = (— l) n -£^F-i(w) is a holomorphic function in u; around w = 0. Involution 
(|6.13p then gives 

Vn(-v) ~ €n(s(t)) = F n (w). 

This completes the proof of the proposition. □ 

As we have noted in Section [5l the residue calculations appearing in the Bouchard-Marino 
recursion formula (|4. 16[) are essentially evaluations of the product of £-forms at the point t 
and s(t) on the Lambert curve, if we truncate the result to the polynomial part. In terms 
of the v-coordinate, these two points correspond to v and —v. Thus we have 
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Corollary 6.3. The residue polynomials of (|5.3p are given by 



(6.15) P a ,b(t)dt 



1 



dt 



ts(t) 
t - s(t) t 2 (t - 1) 



S(t)) + |a+l («(*))&+! (t) 



r] a+1 (v)r] b+1 (v) 
V-i(v) 



vdv 



v=v(t) 



where the reciprocal of 
is defined by 



OO / OC 1 \ 

£ Sfc (2(-l -*)-!)!! = -vll + £(-l)* Sjfc _l 

fc=0 V k=l V ; " / 



\ m =0 U=l v 



.2/c 



)!! 



Proof. Using the formulas established in Proposition 16.21 we compute 



dt 



ts(t) 
t-s(t) t 2 {t - 1) 



£a+i(t) - L+i (*(*)) - 4+i («(*)) 



»7a+l(u)»76+l(u) 



g a+ i(t)+^ a+ i(8(t)) jj+ijt) + 6+i («(*)) 

2 2 

'/a+l(«)»76+l(«) - %W%W| (-v)dv 

vd-u + (const + 0(w)^dv. 

- 0. This completes the proof of (|6.15|) . □ 



2?7-i (w) 

From (|2.15p we see [(const + 0(w))dv\ 

For the terms involving the Bergman kernel, we have the following formula. 
Proposition 6.4. As a polynomial in t and tj, we have the following equality: 



(6.16) P n (t, tj)dt ® dtj = d t 



ts(t) ( i n+1 {t)ds{t) £ n+1 (s(t))di 



t - s(t) \ s(t) - tj 



t - tj 



r \ .. finite 
Vn+liVj) _ J_ fVj_ 



m=0 



2m 



■u=«(t) 
«j=«(*j). 



In £/ie we evaluate the expression at v = v(t) and Vj = v(tj), then expand it as a 
series in | and j-, and finally truncate it as a polynomial in both t and tj. 

Proof. From the formulas for £ n (i) and rj n (v), we know that both expressions have the 
same degree 2n + 2 in t and tj. Since the powers of Vj in the summation Yl^=o (if) m ^ s 
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non- negative, clearly we have 



7?_1 (v) v z ^-^ \ V 
1 iw m=0 



u=u(t) 
«i=«(*i). 



0. 



Thus we can replace the RHS of (|6,16p by 



finite 



ifc+i(«)-**a(«i). ^p\-» ( _ u)dl 



m=0 



ri=i;(t) 
«j=«(*i). 



Since the degree of T] n+ 3^v(tj)) in ij is 2n + 1, the finite sum in m of the above expression 
contributes nothing for m > n + 2. Therefore, 



finite 



2m 



dt- 



L ? 



r)-i(v) v 2 ^— ' V v / 

v y m=0 
/£n+l(*) - £rH-l(*j) . ^ri+l(w) - -Pn+l(%') 



t/=i;(t) 



i — \ w — Wj 



ts(t) in+l{t) - £n+l(tj) 



w — w 



(-dw) 



v=v(t) 
Vj=v(tj), 



t — s(t) W — Wj 

because of (16.81) . We also used the fact that 

finite 

£ 



dw) 



v=v(t) 
Vj=v(tj). 



1 '^ e /f 7 \ 2m 1 dw „, 2 . , 

i vdv = h 0(w™ +2 )dw, 



m=0 



2 W — Wj 



and that is holomorphic along w = Wj. Let us use once again — £ n+ i(i) 



£„+i(s(i)) +2F n+ i(w) and 



n2 — 



dv. We obtain 



L ? 



ts(t) £ n+1 (t) - Zn+xitj) 



-dw) 



v=v(t) 

Vj=v(tj). 



t — s(t) W — Wj 

ts(t) f in+i(t) - £ n +l(tj) dv 



d t - 



ds(t) 



+ dt 



t — s(t) \ —v — Vj ds(t) 

ts(t) ( i n+1 (s(t)) - £ n+1 (s(tj)) dv 



t - S 



V — Vj 



dt 



dt 



v=v(t) 

v j= v l t i). 



v=v{t) 



d t - 



ts(t) ( L+i(t) - L+i(tj) s(t)-tj dv(s{t)) 



t-s(t)\ s{t)-tj v(s(t)) - v(tj) ds(t) 



ds(t) 
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+ d u 



ts(t) ( L+i{s{t)) - £ n+l ( S (tj)) 



t - s(t) 



t - ti 



t — tj dv{t) 
v(t) - v(tj) dt 



dt 



Here we remark that 

ts(t) ( in+l(t) ~ L+l(tj] 



dt. 



t - s(t) 



s(t) - tj 



t - 



t - s(t) \ s(t) - tj 



t - ti 



because the extra terms in the LHS do not contribute to the polynomial part in t. Therefore, 
it suffices to show that 



(6.17) d t 



ts(t) j n+1 (t) - j n+1 (tj) ( s(t) - tj dv(s(t)) 

\v(s{t))-v(tj) ds(t) 



t-s(t) S(t)-tj 



1 ds(t) 



+ <k 



ts(t) &H-1 (*(<)) (*(**)) 



d t 



+ dt 



t - s(t) 

ts(t) 
t - s(t) 
ts(t) 
t - s(t) 

ts(t) 
t - s(t) 

-dt. 



t - ti 



t - ti 



dv{t) 



in+l(t) - £,n+l(tj 
£ n+l (s(t)) - £ n+1 (s(tj)) 

£n+l(*) _ £n+l(tj 
ts(t) 



v(t) - v(tj) dt 
-dv(t) ds(t) 



1 dt 



v(t) - V{tj) s(t) - tj 
dv(t) dt 



v(t) — v(tj) t — tj 

dv(t) ds(t) 



£n+l(*) - Sn+l(tj) 



v(t)-v{tj) S{t)-tj 

dv(t) dt 



t - s(t) 

in light of (|6,8p . At this stage we need the following Lemma: 
Lemma 6.5. For every n > we have the identity 

—dv ds(t) dv 



v(t) — v(tj) t — tj 



0, 



+ 



(6.18) 



= df 



-V - Vj S(t) -tj V- Vj 



+ 



dt 



t - U 



v=v(t) 
«i=«(*i). 



Proof of Lemma. First let us recall that B(t,tj) = dt - y-j—j-j is the Bergman kernel of the 

Lambert curve C, which is a symmetric quadratic form onCxC with second order poles 
along the diagonal t = tj. The function v = v(t) is a local coordinate change, which 
transforms v = to t = oo. Therefore, the form ^ — j^- is a meromorphic 1-form locally 
defined on C x C, which is actually holomorphic on a neighborhood of the diagonal and 
vanishes on the diagonal. Therefore, it has the Taylor series expansion in j and without 
a constant term. 

Since v(s(t)) = -v(t), the form jg. - - ^- - 11 ~ f * 



S{t)-tj 



is the pull-back of — 7^- via the 

local involution s : C —* C that is applied to the first factor. Thus this is again a local 
holomorphic 1-form on C x C and has exactly the same Taylor expansion in and j-. 

Therefore, in the j— expansion of the difference . — f. w . — 77^77- + 7^77, each coefficient 
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does not contain a constant term because it is cancelled by taking the difference. It implies 
that the difference 1-form does not contain any terms without \. In other words, we have 



(6.19) 







-dv 



ds(t) 



dv 



+ 



dt 



-V — Vj s(t) — tj V — Vj t — t 



v=v(t) 
Vj=v(tj). 



Note that we have an expression of the form 

—dv ds(t) dv dt 



(6.20) 



+ 



v=v(t) 



-V — Vj s(t) — tj V — Vj t — tj 

where / is a power series in one variable and F a power series in two variables. Therefore, 

—dv ds(t) dv dt 



(6.21) 



= dt 



f" 



-V — Vj s(t) — tj v — V, 



+ 



t - L 



v=v(t) 

Vj=v(tj). 



Lemma follows from (I6.19P and (I6.21|) . 
It is obvious from ([67T9]) and (I6T20D that 



□ 



(6.22) = d t 



t r 



ts(t) 



-dv 



ds{t) 



dv 



+ 



dt 



3 t — s(t) \ —V — Vj s(t) — tj V — Vj t — t 



v=v(t) 
Vj=v(tj). 



Since £ n +i(i) is a polynomial in t, (16.170 follows from (|6.22[) . This completes the proof of 
the proposition. □ 

7. Proof of the Bouchard-Marino topological recursion formula 

In this section we prove (|1.5|) . Since it is equivalent to Conjecture 14.81 we establish the 
Bouchard-Marino conjecture. Our procedure is to take the direct image of the equation 
(|3.12p on the Lambert curve via the projection ir : C — » C. To compute the direct image, it 
is easier to switch to the coordinate v of the Lambert curve, because of the relation (|2.17p . 
This simple relation tells us that the direct image of a function f(v) on C via the projection 
7r : C — * C is simply the even powers of the v- variable in f(v): 

W = /(«) + /(-«)• 

After taking the direct image, we extract the principal part of the meromorphic function in 
v, which becomes the Bouchard-Marino recursion (|1.5p . To this end, we utilize the formulas 
developed in Section [6l 

Here again let us consider the £ = 1 case first. We start with Proposition 13.41 

Theorem 7.1. We have the following equation 
(7.1) - £ (r n A^l)) g ^ 1 (v)r )n+1 (v) = l £ 



n<3g-2 



a+b<3g-4 



stable 



91+92=9 



r} a+ i(v)r) b+1 (v) +O w (l] 



where O w (l) denotes a holomorphic function in w 
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Proof. We use (j6.8|) to change from the t- variables to the variables. The function factor 
of the LHS of (13.61) becomes 



(2g- l)rj n {v) +ij n+1 (v) - ij^ 1 {v)r] n+1 (v ) +vf L (w) + const + 0(w), 

where fi^w) is a Laurent series in w. The function factor of the RHS is 

ila+i(v)rj b+1 (v) + vf R (w) + const + 0(w), 

where fn(w) is another Laurent series in w. We note that the product of two r/„-functions 
is a Laurent series in w. Therefore, extracting the principal part of the Laurent series in w, 
we obtain 



- £ {T n ^ g (l)) g ^ l (v)Vn + i{v) = \ E 

n<3g-2 a+b<Sg-4 

+ E foWW^ww 

9 1+92=9 

which completes the proof of the theorem. □ 
Corollary 7.2. The cut-and-join equation (I3.5P for the case of £ = 1 implies the topological 



(r a 7iAY_ 1 (l)) 12 



r? a+ i(i;)r/ 6+ i(w) + const + 0(u/ 



recursion 



Proof. Going back to the t-coordinates and using (|2.10p . (I2,ll|) . and (|6.8I) in (17. ip . we 
establish 



(7.2) £ (r„A 3 v (l))^n(t) = i E ((TanA^ 

n<3g-2 a+fe<3g-4 V 

stable \ 



91+92=9 



since 



const + O(to) 
7/-1 (v) 



u=i;(t) 



(1)) 9 -1,2 

?7a + l(f)??fc + l(f) 

0. 



u=u(t) 



From Corollary 16.31 we conclude that (|7.2p is identical to (|5.6p . This completes the proof 
of the topological recursion for 1=1. □ 



We are now ready to give a proof of (jl.5p . The starting point is the Laplace transform of 
the cut-and-join equation, as we have established in Theorem 13.71 Since we are interested 
in the principal part of the formula in the u-coordinate expansion, in what follows we ignore 
all terms that contain any positive powers of one of the Uj's. 

First let us deal with the unstable (0, 2)-terms computed in (|3.9|) . Using (|6.9p . we find 



»7-i(«i) — »7-i(«j) 

modulo holomorphic functions in Wi and Wj . Therefore, the result of the coordinate change 
from the i-coordinates to the ^-coordinates is the following: 



(7.3) E ( T nA V g iX)) g ,l ( (25 - 2 + l)Vn L K) + E Vm+i(vi)vL\{i}(v L \ { 



II L 



1=1 



.'SO 
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-^2v-i{ v i)Vm+i(vi)r]n LW} (vL\{{}) 
i=i J 

1 1 ( 

i=l n L \ii\ a ,b \ 



stable \ 

+ Yl ( r « r ™/ A 9i( 1 )) 9 i,|/|+l^ T ™.7 A 9 2 ( 1 )) 92 ,|j|+i JVa 
9i+g2=g / 



+1 (vi)f]b+i (vi)Vn L ^ } (v L \ {i} ; 



91+92=9 
lUJ=L\{i} 



1 

+ ?YY Y Y^L^^KWgJ-lVnL^yivL^ij}) 

rj m+ i(vi)r]o(vi) - •q m+ x{v j )r]Q{vj) 



2 

t=l j^i n L \ {iJ} 



r)-i{vi) - V-l(vj) 

again modulo terms containing any holomorphic terms in any of w^s. At this stage we 
take the direct image with respect to the projection ir : C — > C applied to the v\ -coordinate 
component, and then restrict the result to its principal part, meaning that we throw away 
any terms that contain non-negative powers of any of the v^s. Thanks to (j6.8l) . only those 
terms containing rj a (vi)rjb(vi) survive. The last term of (|7.3[) requires a separate care. We 
find 

1 ( ^m+l^O^O^l) ~ Vm+l{Vj)Vo(. v j) Vm+l{vi)Vo(vi) ~ V m +l(Vj)m(.Vj) 



2k 



2 V »7-i(wi) - V-lfa) -V-i(vi) ~ V-i( v j) 



r]-i(vi) 2 - r]-i(vjY z 



_ r) m+ i(vj) _ r] m+1 (vj) ST^fvi 
vf-v] vl j^W 

modulo terms containing non-negative terms in vj. Thus by taking the direct image and 
reducing to the principal part, (|7.3p is greatly simplified. We have obtained: 

Theorem 7.3. 



2 Y Y ( (Wolu^s-iWIj-M+i 

n L\{i} a,b \ 



stable \ 

+ Yl (TaTnAl(l)) 9lM+ l(T b TnAUl)) 92 ,\J\ + l 
31+92=9 / 



Va+l(vi)Vb+l(vi)Vn LX{1} (v L \ { i } ) 



IUJ=L\{1} 

, « finite , -,2k 

+ \Y Y E^^^a^i))^^ 

i>l«L\{i,j} m 1 k=0 v 17 

modulo terms with holomorphic factors in v^. 
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We note that only finitely many terms of the expansion contributes in the last term of 
(|7.4p . Appealing to Corollary 16. 31 and Proposition 16,41 we obtain (II. 5ft . after switching back 
to the t-coordinates. We have thus completed to prove the Bouchard-Marino conjecture [3]. 

Appendix. Examples of linear Hodge integrals and Hurwitz numbers 



In this Appendix we give a few examples of linear Hodge integrals and Hurwitz numbers 
computed by Michael Reinhard. 



g = 2 


£=l 


(73^1)2,1 


1 

480 






1 = 2 


<t|Ai) 2 ,2 


5 
576 






g = 3 


1 = 1 


(76Al} 3) l 


7 


(75^2)3,1 


41 




138 240 


580 608 


1 = 2 


(•7275^1)3,2 


323 


(t 2 T 4 A2)3,2 


2329 




483 840 


2 903 040 


(7374^1)3,2 


19 


(7"3 2A 2>3,2 


1501 




17 920 


1451520 


£ = 3 


(T2 r 4^l)3,3 


541 


(t27-3 3 Ai) 3) 3 


89 
7680 


(TfT 3 A 2 )3,3 


859 


60 480 


96 768 


£ = 4 


(t 2 3 T3Ai) 3 ,4 


395 
3456 


(t|A 2 ) 3 ,4 


17 
192 




g = A 


l=\ 


(^9^1)4,1 


1 


(^2)4,1 


1357 


(7^3)4,1 


13 


1 244 160 


696 729 600 


6 220 800 


£ = 2 


(r 2 r 8 Ai)4 i 2 


841 


(t2T 7 A 2 )4,2 


33 391 


(7"37"5A 3 )4,2 


2609 


38 707 200 


696 729 600 


29 030 400 


(73 r 7 Ai) 4i 2 


221 


(t"3T6A 2 )4,2 


1153 


( r 4 2 A3>4,2 


6421 


4147 200 


11059 200 


58 060 800 


(747^1)4,2 


517 


(t"4T5A 2 )4,2 


979 




5 806 080 


6 451200 


(7-5^1)4,2 


1223 


(72^3)4,2 


5477 




11612160 


116121600 


£ = 3 


(r 2 2 r 7 Ai) 4 ,3 


3487 


(t3T4Ai} 4 ,3 


137 


(rf 7^2)4,3 


58 951 


5 806 080 


46 080 


16 588 800 


(T2T3r 6 Ai)4 )3 


50 243 


(T|r 6 A 2 )4,3 


137 843 


(f2 2 7^3)4,3 


241 


38 707 200 


116121600 


230 400 


(t2T4T 5 Ai)4 i 3 


2597 


(t2T 3 T5A2}4,3 


577 


(72737^3)4,3 


27 821 


1 382 400 


258 048 


16 588 800 


(r|r 5 Ai)4,3 


3359 


(t 2 t| A 2 )4,3 


2657 


(T3 3 A3>4,3 


4531 


1 382 400 


967 680 


2 073 600 


g = 5 


£= 1 


(Ti2Al)5,l 


1 


(7^3)5,1 


71 




106168 320 


1114 767 360 


(tuA 2 )5,1 


577 


(t"9A 4 )5,1 


21481 




16 721510 400 


367 873 228 800 



Table 1. Examples of linear Hodge integrals. 



Some examples of g = 5 Hurwitz numbers: 

^5,(1) =0 h 5> (4) = 272 097 280 
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h 5{2) = 1/2 /i 5(5) = 333 251 953 125 

/j 5(3) = 59 049 h 5i{6) = 202 252 053 177 720 



1 


9=1 


g = 2 


g = 3 


g = 4 


(1) 














(2) 


1/2 


1/2 


1/2 


1/2 


(1,1) 


1/2 


1/2 


1/2 


1/2 


(3) 


9 


81 


729 


6561 


(2,1) 


40 


364 


3280 


29 524 


(1,1,1) 


40 


364 


3280 


29 524 


(4) 


160 


5824 


209 920 


7 558 144 


(3,1) 


1215 


45 927 


1673 055 


60 407127 


(2,2) 


480 


17472 


629 760 


22 674 432 


(2,1,1) 


5460 


206 640 


7 528 620 


271 831560 


(1,1,1,1) 


5460 


206 640 


7 528 620 




(5) 


3125 


328125 


33 203 125 


3 330 078 125 


(4,1) 


35 840 


3 956 736 


409 108 480 


41394 569 216 


(3,2) 


26460 


2 748 816 


277118 820 


27 762 350 616 


(3,1,1) 


234 360 


26184060 


2 719 617120 


275 661886 500 


(2,2,1) 


188160 


20160 000 


2 059 960 320 


207 505 858 560 


(2,1,1,1) 


1 189440 


131670 000 


13 626 893 280 




(1,1,1,1,1) 


1 189440 


131670 000 






(6) 


68 040 


16 901 136 


3 931876 080 


895132 294056 


(5,1) 


1 093 750 


287109 375 


68 750 000 000 


15 885 009 765 625 


(4,2) 


788480 


192 783 360 


44490 434 560 


10 093 234 511360 


(4,1,1) 


9 838 080 


2 638 056 960 


638 265 788 160 


148 222 087453440 


(3,3) 


357210 


86113125 


19 797948 720 


4487187 539 835 


(3,2,1) 


14 696 640 


3 710 765 520 


872 470 478 880 


199 914163 328 880 


(3,1,1,1) 


65 998 800 


17634 743 280 


4 259 736 280 800 




(2,2,2) 


2 016 000 


486 541 440 


111644332 800 


25 269 270 586 560 


(2,2,1,1) 


80438400 


20 589 085 440 


4874 762 692 800 




(2,1,1,1,1) 


382 536 000 


100 557 737 280 






(1,1,1,1,1,1) 


382 536 000 









Table 2. Examples of Hurwitz numbers for 1 < g < 4 and \/i\ < 6. 
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